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We observe a spontaneous parity breaking bifurcation to a ferromagnetic state in a spatially- 
trapped exciton-polariton condensate. At a critical bifurcation density under nonresonant excita¬ 
tion, the whole condensate spontaneously magnetizes and randomly adopts one of two elliptically- 
polarized (up to 95% circularly-polarized) states with opposite handedness of polarization. The 
magnetized condensate remains stable for many seconds at 5 K, but at higher temperatures it can 
flip from one magnetic orientation to another. We optically address these states and demonstrate 
the inversion of the magnetic state by resonantly injecting 100-fold weaker pulses of opposite spin. 
Theoretically, these phenomena can be well described as spontaneous symmetry breaking of the spin 
degree of freedom induced by different loss rates of the linear polarizations. 


Condensation of exciton-polaritons (polaritons) spon¬ 
taneously breaks the global phase symmetry [1-5]. Ow¬ 
ing to their easy optical interrogation, high-speed (ps) 
interactions, and macroscopic coherence (over hundreds 
of microns) [6], polariton condensates are excellent candi¬ 
dates to probe and exploit for sensing [7, 8], spinoptron- 
ics [9-11], new optoelectronic devices [12-14], and quan¬ 
tum simulators [15]. The driven-dissipative multicom¬ 
ponent polariton system can undergo additional bifurca¬ 
tions and condense into states which are not eigenstates 
of the single-particle Hamiltonian, but many-body states 
with reduced symmetry [16, 17]. Thus, we should ex¬ 
pect that two-component exciton-polariton condensates 
can also show spontaneous symmetry breaking bifurca¬ 
tions in their polarization state. Spin studies of micro¬ 
cavity polaritons have been of great interest in recent 
years [18-29]. However, spontaneous symmetry-breaking 
bifurcation of spin has not been observed before. 

Here, we demonstrate spontaneous magnetization in an 
exciton-polariton condensate, as a direct result of bifurca¬ 
tions in the spin degree of freedom. Utilizing an optically 
trapped geometry, condensates spontaneously emerge in 
either of two discrete spin-polarized states that are sta¬ 
ble for many seconds, > 10^^ longer than their formation 
time. These states emit highly circularly-polarized co¬ 
herent light (up to 95%) and have opposite circular po¬ 
larizations. The condensate stochastically condenses in 
a left- or right-circularly polarized state, with an occur¬ 
rence likelihood that can be controlled by the ellipticity 
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of the nonresonant pump. The two spin-polarized states 
can be initialized and switched from one state to another 
with weak resonant optical pulses. Our system has po¬ 
tential applications in sensing, optical spin memories and 
spin switches, and it can be implemented for studying 
long-range spin interactions in polariton condensate lat¬ 
tices. 

This article is structured as follows: in Section I we re¬ 
view trapped polariton condensates and the current un¬ 
derstanding of polarization in untrapped polariton con¬ 
densates. In Section H we present the key theme of this 
work, which is the spontaneous buildup of stochastic cir¬ 
cular polarization. In Section HI we propose a theoretical 
framework for the phenomena discussed in this work. We 
show that stochastic circular polarization is a signature 
of spontaneous parity breaking. In Section IV we time- 
resolve the coherent driving of the spin, with resonant 
excitation. We furthermore investigate the stability of 
the spin-polarized states against thermal noise and con¬ 
clude in Section V. 


I. INTRODUCTION 

Exciton-polaritons are spinor particles formed by the 
strong coupling of excitons in a semiconductor quan¬ 
tum well with photons in the microcavity in which they 
are embedded [30]. We create optically-trapped polari¬ 
ton condensates by nonresonant excitation of a semicon¬ 
ductor microcavity membrane [see FIG. 1(a) and Ap¬ 
pendix A]. The excitation beam is shaped into a 4- 
spot pattern [shown by dashed circles in FIG. 1(b)]. 
The short-wavelength continuous wave (GW) linearly- 
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FIG. 1. (a) Scanning electron microscopy image of a mem¬ 
brane microcavity, cleaved across the middle. The membrane 
allows resonant excitation from the back side of the cavity, 
(b) Above threshold hot polaritons condense into a single co¬ 
herent state located at the center of the trap (pump spots 
shown by dashed circles), which is found at k\\ =0 (c). The 
spin state of the condensate is determined by studying the 
emission polarization. 


polarized pump injects an electron-hole plasma at each 
pump spot, which rapidly relaxes to form excitons, in 
the process losing all phase information. These reser¬ 
voir excitons then scatter into polariton states via multi¬ 
ple phonon-polariton and stimulated polariton-polariton 
collisions [31], and they feed the zero-momentum ground 
state at the center of the trap. Because of their large 
effective mass, excitons typically diffuse only very small 
distances and stay within 1 pm of the pump spots. Micro¬ 
cavity polaritons, however, are 10,000 times lighter giv¬ 
ing longer diffusion lengths. Driven by their repulsive 
excitonic interactions, polaritons can thus travel large 
distances away from the pump spots within their life¬ 
time [32]. Once the density inside the trap exceeds the 
condensation threshold, a macroscopically coherent con¬ 
densate is formed [FIG. 1(c)]. Because the condensate 
overlaps only weakly with the pump spots, it shows a 
narrower linewidth and less decoherence than unconfined 
condensates [33]. The optical trapping method used here 
is similar to optical lattices in cold atomic systems [34], 
but with the major difference that the optical potential 
also provides gain [32, 35, 36]. 

Polaritons in quantum-well microcavities have two 
Jz = (spin-up or spin-down) projections of their total 
angular momentum along the growth axis of the struc¬ 
ture, which correspond to right- and left-circularly polar¬ 
ized photons emitted by the cavity, respectively. When 
the excitation is linearly-polarized an equal population 
of spin-up and spin-down excitons forms in the reservoir. 
An initially spin-balanced reservoir, in the absence of pin¬ 
ning to any crystallographic axis, is expected to give a 
condensate with a stochastic linear polarization [37]. In 


most experiments with polariton lasers the condensates 
have been found to be linearly polarized along one of the 
crystallographic axes [1, 2]. Nevertheless, in some cases 
a circularly polarized polariton lasing has also been ob¬ 
served [5, 19, 23, 38, 39]. Formation of a circularly polar¬ 
ized condensate is usually associated with the effects of 
TE-TM splitting, or bosonic amplification of the seed po¬ 
larization of condensates. In all these cases, a circularly 
polarized condensate is observed when the symmetry be¬ 
tween the spin-up and spin-down polaritons has been ex¬ 
plicitly broken in some fashion, by either the pumping 
geometry (seeding with a circularly polarized pump) or 
by the imposed rotation of the Stokes parameters due to 
polarization splitting. As a result, the observed circu¬ 
lar polarization was never stochastic (i.e. it is fixed each 
time the condensate is excited, but it is different on each 
realization). 


II. SPONTANEOUS BUILDUP OF CIRCULAR 
POLARIZATION 


Our experiments reveal a completely different behav¬ 
ior to the previous reports on the polarization of polari¬ 
ton condensates. We observe a strong degree of circu¬ 
lar polarization 60% < \sz\ < 95% [FIG. 2(a)] when 
the excitation is linearly polarized (to better than 1 in 
10^), which is stable for many seconds [FIG. 2(b)]. The 
condensate stochastically adopts either of two opposite 
circular-polarization states in each realization of the ex¬ 
periment. We call these two states the spin-up (5>^) and 
spin-down (s^) states. By mapping the polarization of 
the photoluminescence (PL) we measure the polarization 
vector Sx^y^z = 

I is the measured intensity for horizontal (iL), vertical 
(U), diagonal (D), anti-diagonal (A), right-circular (O) 
and left-circular (O) polarizations. We measure all com¬ 
ponents of the polarization vector (pseudospin) simulta¬ 
neously and plot the mean of the spin-up and spin down 
states on the Poincare sphere separately for 1000 realiza¬ 
tions [FIG. 2(c)]. In each realization the wavefunction 
of the condensate spontaneously collapses into one of the 
two discrete spin-polarized states which have opposing 
circular and diagonal components, marked by blue and 
orange vectors in FIG. 2(c). The linear axis along which 
the pseudospin flips (marked here as diagonal) does not 
depend on the geometry of the trap, and it changes direc¬ 
tion with the position of the condensate on the sample. 

To demonstrate that the buildup of circular polariza¬ 
tion is truly spontaneous, we illuminate the sample with 
a long-duration pulse of 9 ms using an acousto-optic mod¬ 
ulator (AOM). A condensate builds up and picks a ran¬ 
dom state (e.g. s^) and stays in that particular state as 
long as the pump pulse lasts [FIG. 3(a)]. We then repeat 
the same measurement but this time we modulate the 
pump intensity so that a condensate is created and de¬ 
stroyed every 2ps. FIG. 3(b) shows how the condensate, 
although being stable for many seconds, picks a random 
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FIG. 2. (a) Polarization-resolved spatial image of the only two realizations observed in a single trapped condensate in a 4-spot 

trapping geometry, (b) At T = 5K the spin-polarized states remain stable for many seconds, (c) Simultaneously measured 
components of the pseudospin for 1000 realizations. The total degree of polarization is 0.93 ± 0.03 and the two average 
pseudospin states are = [-0.22,0.19,-0.94] (blue vector) and Sf = [-0.22,-0.14,0.96] (orange vector). The measurement 
error for each component is < 5%, and the variance is ~ 1%. 


polarization at every successive realization due to ran¬ 
dom initial conditions at the onset of condensation. In 
this case the parity symmetry is broken spontaneously. 
The buildup of circular polarization is independent of the 
precise position on the sample, of the sample orientation, 
of the pump spot orientation, and of the polariton detun¬ 
ing. 

We can also explicitly break the symmetry in our ex¬ 
periments by changing the ellipticity of the pump laser 
and measure the contrast of the occurrence frequency of 
spin-up (/|) to spin-down (/;) condensate realizations, 
^ = (/t - h)/{h + /t)- The probability of a realiza¬ 
tion resulting in state is equal to (1 ± (^)/2 respec¬ 
tively. FIG. 3(c) shows ^ as a function of the pump 
circular polarization (sz^pump) averaged over 1000 real¬ 
izations for each point. As the ellipticity of the pump 
is increased from linear {sz^pump = 0) to right circular 
{sz,pump > 0) the probability of creating a condensate 
in the spin-up state increases; conversely, for the left- 
circularly polarized pump {sz,pump < 0), the probability 
of creating a spin-down condensate increases. With a 
linearly-polarized pump we have an equal probability of 
creating a spin-up or spin-down condensate. Although 
the pump laser is nonresonant with the final polariton 
states, the initially created carrier spin is not entirely ran¬ 
domized during their multiple carrier-carrier and exciton- 
phonon scatterings [5, 28, 39, 40]. As a result of this in¬ 
complete spin relaxation of the excited carriers, changing 
the ellipticity of the pump breaks the symmetry of the 
condensate toward the same circular polarization as that 
of the pump. For pump circular polarizations far greater 
than 0.1 the condensate is formed deterministically in the 
same polarization state as that of the pump [FIG. 3(d)]. 

An interesting question here is why we observe the 
spontaneous buildup of circular polarisation, as opposed 
to the linear polarisation that is widely reported in the 
literature [1, 2]. The key difference between the exper¬ 
iments presented in this work and other studies of po¬ 


lariton condensation lies in the excitation geometry. For 
our trapped condensates, the pump and condensate are 
spatially separated, which critically reduces the contam¬ 
inating interactions between the condensate and reser¬ 
voir. The large interaction between untrapped conden¬ 
sates and the unpolarized exciton reservoir results in 
spin-flip scattering of polaritons with reservoir excitons. 
If there is a depolarized reservoir on top of the conden¬ 
sate, the spin-flip scattering processes minimize any im¬ 
balance between circular components of the condensate. 
This minimization leads to quenching of the buildup of 
circular polarization, forcing the polaritons to condense 
only with linear polarization. Moreover it has been shown 
previously that, because of a smaller overlap with the 
reservoir, trapped condensates have a smaller linewidth 
than untrapped condensates [33]. Our careful studies 
with temperature and our theoretical calculations (see 
Section IV D), show how spin noise in the system re¬ 
sults in spin flipping of the condensate. The spin flip 
rate scales exponentially with noise. Larger linewidth 
untrapped condensates have higher spin flip rates, which 
wash out circular polarization effects observed here. 

To show the crucial role of the reservoir excitons we 
place a weak nonresonant linearly-polarized probe beam 
on top of the condensate as shown in FIG. 4(a). The 
probe beam which has just a small fraction of the 4-spot 
pump power (< 0.025 Pth)i induces a reservoir of exci¬ 
tons that overlap with the condensate but it does not 
stop the condensation or reduce the condensate density 
below the critical circular polarization density. We then 
measure the absolute average circular polarization de¬ 
gree in each realization, as a function of the weak probe 
power. As shown in FIG. 4(b), the circular polariza¬ 
tion degree decreases monotonously as the probe power 
increases, demonstrating the quenching of circular polar¬ 
ization due to the influence of the overlapping reservoir. 
This measurement evidences the importance of the sep¬ 
aration of the pump-induced reservoir and the trapped 
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FIG. 3. (a) Emitted circular cr+ and g- intensities for a 

condensate realization when P > Pc , created by a long pulse 
(marked by P). The condensate randomly picks the state 
here. Note the opposite axis directions for cr+ and a- com¬ 
ponents. (b) Same as (a), with pump intensity modulated by 
square wave (shown on top by P). A cut from the middle of a 
9 ms exposure is shown, where each square pulse corresponds 
to a single realization. The condensate randomly picks or 
Si states (marked by grey rectangles), (c) Occurrence con¬ 
trast ^ of spin-up and spin-down realizations measured in (b) 
vs the ellipticity of the pump, (d) Measured circular polar¬ 
ization of condensate vs pump ellipticity. Dashed lines mark 
the region where the condensate initializes stochastically in 
spin-up or spin-down states in each realization. 


(a) 


pump beams 

/\ 




weak 

nonresonant 

probe 



Weak probe power 


FIG. 4. (a) Nonresonant linearly polarized weak probe over¬ 
lapped with the condensate, (b) Average circular polarization 
degree vs power of the probe beam. The blue line is a guide 
to the eye. 


condensate in the observation of ferromagnetic conden¬ 
sates. 

Four scenarios might explain the buildup of the 
stochastic circular polarization depending on populations 
of same- and opposite- (or cross-) spin polaritons: (i) a 
higher interaction energy for cross-spin than for same- 
spin polaritons in a condensate which is in thermal equi¬ 
librium, (ii) a higher gain from the cross-spin reser¬ 
voir, (iii) density-dependent losses enhanced by cross¬ 
spin condensed polaritons, and (iv) linear polarization 
energy splitting accompanied by a dissipation rate split¬ 
ting which destabilizes the linearly-polarized condensate 
(‘spin bifurcation’). As we will show now, our ex¬ 
perimental data strongly suggest that scenario (iv) is 
the correct explanation. In the first scenario (‘energy- 
minimization’), the condensate free energy is minimized 
when it acquires circular polarization. This situation 
could happen in a condensate in thermal equilibrium if 
the interaction energy of cross-spin polaritons is stronger 
than that of same-spin polaritons, meaning the coexis¬ 
tence of cross-polarized polaritons is not favored [41]. 
However, the interaction with opposite spin polaritons 
is well known to be weaker than, and opposite to, that of 
same-spin polaritons [29, 42, 43]. Moreover, being exter¬ 
nally driven and having short lifetimes, polariton conden¬ 
sates are generally far from thermal equilibrium [44]. We 
note that “energetical” mechanisms cannot explain why 
states with elliptical polarization (i.e. not fully circu¬ 
lar) are formed. In the second scenario (‘cross-gain’), in 
order to acquire a circular degree of polarization, the con¬ 
densate must experience larger gain from the cross-spin 
reservoir than from the same-spin reservoir. However, 
the scattering rate from the cross-spin reservoir into the 
condensate is measured to be significantly smaller than 
the same-spin [18]. The third scenario (‘cross-loss’) re¬ 
quires the condensate loss rate to increase when oppo¬ 
site spins are present. If polaritons are directly excited 
by light, observations have suggested biexciton formation 
can produce such enhanced losses [26]. However, this is 
only significant when the relative spectral detuning of 
cavity and exciton is small (< 2meV), which is far from 
the case here. We observe the spontaneous buildup of cir¬ 
cular polarization throughout the entire +2 to —lOmeV 
detuning range. In the remaining scenario (‘spin bifurca¬ 
tion’), which we present in Section HI for the first time, 
the buildup of circular polarization is caused by small dif¬ 
ferences of the energy and dissipation rates of two orthog¬ 
onal linearly-polarized polariton modes, which is present 
even at /c = 0. 

Strain-induced splitting of the linear components of 
polariton condensates has been demonstrated previ¬ 
ously [45, 46]. In our sample this splitting varies de¬ 
pending on the position on the sample. We observe a 
linear polarization energy splitting of up to ^100 peV de¬ 
pending on the position of excitation on the membranes. 
Note that, all the effects reported here are also seen on 
unetched samples, so strain from patterning is not cru¬ 
cial. However, we observe a higher splitting at the edges 
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of the membrane than in the middle, as there is more 
stress in the structure at the edges. Due to the curvature 
of the cavity stopband, any energy splitting is accompa¬ 
nied by a difference of the linewidth (dissipation rate), as 
explained in Appendix B. This energy splitting between 
the two linear components combined with a difference in 
dissipation rates causes the polarization of the conden¬ 
sate to change from linear polarization to circular at a 
critical density (see Section III). 

We emphasize that the stochastic circular polarization 
here cannot be explained in the framework of the opti¬ 
cal spin Hall effect [28, 47]. In our trapping geometry 
the condensate is formed at the ground state with (k = 
0 , Sk = 0.4 pm“^) [FIG. 1 (c)], where the transverse- 
electric and transverse-magnetic (TE-TM) splitting van¬ 
ishes [38]. The trap diameter here is 6 times smaller than 
the observed spin-ring patterns measured for a nonequi¬ 
librium condensate formed at much higher in-plane wave- 
vectors in the same sample [28]. Moreover, the geometry 
or the orientation of the trap does not affect the polariza¬ 
tion state of the condensate. We also see the stochastic 
circular polarization with a ring-shaped trap and also in 
high-order spatial mode condensates [48, 49]. Finally, it 
should be noted that any theoretical picture that assumes 
the buildup of circular polarization arises only because of 
the geometrical arrangement of the pump, would neces¬ 
sarily fail to explain the most essential part of this work, 
which is the spontaneous symmetry breaking (stochastic 
behavior). 


III. SPIN BIFURCATION THEORY (BROKEN 
PARITY) 


Our theory is a development of the theory of polariton 
weak lasing in two coupled condensation centers [16] now 
for the case of the spin degree of freedom. Here, we 
have right and left circular polarizations instead of two 
separated condensates, and we also allow for the gain- 
saturation nonlinearity in the system. 

The order parameter for an exciton-polariton conden¬ 
sate is a two-component complex vector T = [t/^+i, ' 0 _i]^, 
where and are the spin-up and spin-down wave 
functions. The components of the order parameter define 
the measurable condensate pseudospin S = ( 1 / 2 ) (T’*’ • 
tr • T), and the normalized spin vector s = S/S, where 
^x,y,z are the Pauli matrices. The components of this 
vector contain information about the intensities and rel¬ 
ative phases of the emitted light. The order parameter 
evolves according to the driven dissipative equation 

- 5 (7 - 

+ I [(ai + a2)S + (tti — a2)SzO'z] 


or in components: 

V’+I = - yiS)^+i - |(7 - i£^)V’-i 

- i(ai|V’+iP +a2|V’-iP)^+i, 
'<P-i = - y{S)'ip-i - i(7 - ie)'ip+i 

- + a 2 |t/'+ip)^/'-i. 


Here, g{S) = T — W -\-rjS is the pumping-dissipation bal¬ 
ance, F is the (average) dissipation rate, W is the incoher¬ 
ent in-scattering (or ‘harvest’ rate), and r] captures the 
gain-saturation term with S = (IV^+iP + |'0_ip)/2 [50]. 
Here, this gain saturation depends on the total occupa¬ 
tion of the condensate (treated more generally in Ap¬ 
pendix C). It is assumed now that X (horizontal) and Y 
(vertical) linearly-polarized single-polariton states have 
different energies and dissipation rates. The energy of 
the X-polarized state is shifted by — 5 / 2 , and the en¬ 
ergy of the Y-polarized state by -h 5 / 2 . The dissipation 
rate from the X-polarized state is F + 7 , while the dis¬ 
sipation rate from the Y-polarized state is F — 7 (see 
also Appendix B). Finally, ai is the repulsive interaction 
constant for polaritons with the same spin, and 0^2 is the 
interaction constant for polaritons with opposite spins. 

From Eq. 1 we obtain for the components of the pseu¬ 
dospin vector (a = 0^1 — 0 ^ 2 ): 


■g{S)S,-jS-aSzSy, 

(3a) 

+ OtSzSx, 

(3b) 

■g{S)Sz - eSy, 

(3c) 


and the related equation for the total spin S = —g{S)S — 
There are two sets of solutions, which we call here 
the paramagnetic and ferromagnetic solutions. 

a. Paramagnetic solutions. These give simple con¬ 
densation into either X or Y linearly-polarized states. 
The Y state possesses the longest lifetime, and the con¬ 
densation threshold is reached for this state first at 
Wi = F — 7 . There is no parity breaking for this con¬ 
densate: Sy = Sz = 0 , Sx = —S with S = {W — Wi)/r], 
so that the occupations of -hi and —1 components are 
equal. However, this condensate solution becomes unsta¬ 
ble for W > W 2 . The values of the critical occupation Sc 
and the critical pumping rate are 

2 2 

Se = W2 = Wi+ vSc. (4) 

as 

Note that this instability is present also for equal dissi¬ 
pation rates, i.e. when 7 = 0 . In this case, the system 
(Eq. 3a-c) describes the self-induced Larmor precession of 
the pseudospin vector. Incorporating energy relaxation 
(e.g., using small negative 7 ) then leads to the formation 
of the X-polarized condensate—an intuitively expected 
result. 

b. Ferromagnetic solutions. The key ingredient of 
our theory is the presence of the 7 > 0 parameter de¬ 
scribing the variation of dissipation rates. This parame¬ 
ter allows the formation of the ‘weak lasing’ regime [16], 
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FIG. 5. (a) Numerical calculations using Eq. 3 for the case when the condensate occupation is below the critical occupation 

Sc (dotted lines) and when above Sc (solid lines). Blue lines show the occupation of the spin-down component and orange 
lines show the spin-up component, (b) Numerical calculations of the Stokes parameters at different powers using Eq. 3. (c) 
Measured power dependence of the condensate Stokes parameters, (d) 2D simulations, with asymmetric initial conditions 
[sz{t = 0) < 0.01] show density of the components of the wave function. Pump spots are marked by dashed lines (see 
Appendix G for parameters). 


which is characterized by two important features: (i) the 
X-polarized condensate is also unstable, and (ii) when 
the F-polarized condensate loses stability at the critical 
occupation S'c, it continuously transforms into one of the 
two ferromagnetic states. While Eqs. (2) are parity sym¬ 
metric, i.e., they are not affected by the interchange of 
left and right circular polarization, the new solutions are 
characterized by broken parity symmetry and by sponta¬ 
neous formation of either left or right elliptical polariza¬ 
tion. These solutions are 


. _ 9iS)^ 


Sy = - — S^, (5a) 


5 . 



7 ^ 

+9{S)‘^ 


l[e^+9{Sr] 

asg{S) 


(5b) 


where the positive root of the second equation in Eq. (5b) 
should be taken. We note that while the sign of Sx is al¬ 
ways negative, Sy and Sz have opposite signs for the two 
solutions. This means that the left-circular component 
is accompanied by a diagonal component, and the right- 
circular by anti-diagonal. Moreover, if these components 
change for some reason, they mirror each other as long 
as the total condensate occupation stays fixed. We label 
these two solutions as the and spin states. 

The spin-independent model for the gain saturation 
used in this Section is sufficient to describe the experi¬ 
mentally observed features. However, the spin relaxation 
in the reservoir can be slow [40], and in this case the 
model should be modified to allow the saturation terms 
to depend on the individual occupations of the left- and 
the right-circular polarization components rather than 
the total occupation only. The parity breaking is still 
present after this modification; however, the stability of 
solutions becomes more complex. The ferromagnetic so¬ 
lutions can now become unstable and transform into pe¬ 


riodic cycles [51], the dynamics of the pseudospin can be¬ 
come irregular, and this can also result in the formation 
of the stable X-polarized condensate at high pumping 
powers. See Appendix C for more details. 

Numerical calculations for the occupation of the two 
circular components of the wavefunction when S < Sc 
(dotted lines) and when S > Sc are shown in FIG. 5(a). 
Here, the condensate is initialized with a small asymme¬ 
try in spin-up and spin down occupations (<1%). Be¬ 
low the critical occupation Sc^ the condensate is linearly 
polarized, but when the occupation is increased above 
the threshold Sc, the condensate adopts one of two el- 
liptically polarized configurations depending on the ini¬ 
tial conditions. In the experiment, the stochastic be¬ 
havior is due to random spin fluctuations at the onset 
of the condensation. In theory, we can reproduce it by 
randomly setting the initial conditions. Numerical cal¬ 
culations of the condensate polarization versus excita¬ 
tion power are shown in FIG. 5(b). Directly at the con¬ 
densation threshold the condensate is linearly polarized, 
but once it reaches the critical occupation (Pc = 1.3Pth, 
marked by a dashed grey line), the linear component is 
quenched and circular polarization builds up. This be¬ 
havior reproduces the experimental data, as shown in 
FIG. 5(c). We observe an initial buildup and subsequent 
quenching of linear polarization with the continuing in¬ 
crease of circular polarization at Pc = 1.2bPth (marked 
by a dashed grey line, with total intensity marked by a 
dotted red line). Once circular polarization is achieved, 
the orientation of the condensate circular polarization be¬ 
comes stochastic under linearly polarized pumping. 

We can extend Eq. (1) and account for 2D real-space 
degrees of freedom by using complex Ginzburg-Landau- 
type equations [8, 35, 50], which in addition to the pump 
and decay also incorporate a repulsive potential due to 
the excitons in the pump spots and an energy relax- 
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ation [52] for polaritons in the trap: 

+ (1 — iA)|— [(o^i + 0 ^ 2 )^ + (<^1 ~ C^2)Sz^z] ^ 

1 1 

- ^ k (6) 

where m* is the effective mass of the polaritons. The 
harvest rate is given by VF = rP/Tji, where P is the 
spin-independent spatial profile of the excitation, T^^ is 
the decay rate of the exciton reservoir, and r is the in¬ 
coming rate of polaritons into the condensate. The gain 
saturation is given by 7 ^ = r^P/T^. The repulsive poten¬ 
tial due to the interaction of polaritons with the exciton 
reservoir is given by Vp = ^grN P ^QpP, where and 
gp are the interaction constants of polaritons with the 
exciton reservoir and the pump spot respectively, and 
N = g{S)/r is the density of the exciton reservoir. Here, 
A ^ 1 is a phenomenological constant that gives the 
energy relaxation. The density profile of the two circu¬ 
lar components of the wave function in the steady state 
[FIG. 5(d)], for the case of a trapped condensate in the 
middle of the 4 pump spots, exhibits a circular polariza¬ 
tion degree of \sz\ = 0.69. Note that with 7 = 0 and 
only polarization splitting (including TE-TM splitting), 
our 2D simulations do not show bistable condensation. 

It is important to note the differences between the fer¬ 
romagnetic states we discuss here and the magnetiza¬ 
tion transition in equilibrium cold atom systems [53, 54]. 
First, the parity breaking bifurcation described above 
does not reduce the energy of the system (unlike for 
atoms). In fact the energy of elliptically polarized states 
is higher than that of linearly polarized states. Second, 
the in-plane components of the spin do not vanish com¬ 
pletely. Third, we have a magnetized condensate with 
only two possible orientations, whereas in atomic systems 
ferromagnetic domains with continuous variable orienta¬ 
tion are observed. 

If the Hamiltonian and the initial state of a system 
are symmetric under the exchange of spin-up and spin- 
down components, but the final state is not, the parity 
symmetry is spontaneously broken. This is indeed the 
case here: We excite an equal population of spin-up and 
spin-down polaritons, which spontaneously condense, but 
form highly circularly polarized macroscopic states in the 
absence of any external magnetic field. 

IV. RESONANT EXCITATION 

The GaAs substrate commonly used in the fabrica¬ 
tion of GaAs microcavities is opaque at the emission 
wavelength (^ 800 nm) of the cavity polaritons. As a 
result, the back side of the cavity is resonantly inaccessi¬ 
ble. Resonant excitation of the cavity from the front side 
has complications with backscatter from the laser, espe¬ 
cially in high finesse cavities and at normal incidence. 


where the condensate emission mode is located. To cir¬ 
cumvent this problem we chemically etch the substrate 
to form membranes of 8.8 pm thickness and 300 pm di¬ 
ameter [see FIG. 1(a)]. For resonant excitation, we use 
a narrow-linewidth (< 2 GHz) GW laser, which is am¬ 
plitude modulated with a second AOM. We call this 
resonant laser the ‘gate’. We use two photomultipliers 
and a fast oscilloscope to time resolve the left- (cr_) and 
right-circular (cr+) polarization intensity of the conden¬ 
sate emission. The resonant excitation laser, the nonres¬ 
onant pump laser, the cameras and the oscilloscope are 
all synchronized, which allows us to vary the delay time 
and amplitudes of each laser pulse on demand. 


A. Resonant initialization of spin states 



FIG. 6. Measured circular polarization of condensate vs 
pump ellipticity, for three different gate polarizations. Dashed 
lines mark the region where the condensate initializes stochas¬ 
tically in spin-up or spin-down states in each realization. 

We can control the polarization of the trapped con¬ 
densate with the gate laser. In this case, we additionally 
resonantly excite the condensate (which is generated by 
the four nonresonant pump spots) with a second laser 
from the backside of the microcavity membrane. This 
gate can be linearly, left-, or right-circularly polarized. 
FIG. 6 shows the condensate circular degree of polariza¬ 
tion versus that of the pump. In the green curve which 
shows the behavior when the gate is linearly polarized, 
we reproduce the same result as that in FIG 3(d). How¬ 
ever, when a right (or left) circularly polarized gate is 
applied, the curve shifts to the left (or right). In other 
words, with a linearly polarized pump, the condensate 
initializes in a right (or left) circularly polarized state. 
The imbalance caused by the resonant gate cancels out 
with an opposite circularly polarized pump with a circu¬ 
lar polarization degree of |5;s^pump| ^0.1. 


B. Elliptically-polarized pump: coherent driving 

In a first experiment, we initialize a condensate in the 
spin-down state (s^) by making the pump laser slightly 
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FIG. 7. (a) Measured intensity of circular polarization components (cr±) for adding gate pulse G (intensity shown on above 

inset). The condensate initializes in the si state, due to the asymmetric pumping condition s^^pump = —0.1. A gate pulse 
with an opposite circular polarization (cr+) is applied, (b) Numerical calculations for (a), showing that during the gate pulse 
G = the spin-up and spin-down densities oscillate with a tt phase shift, (c) Time evolution of the pseudospin, showing 
how the pseudospin orbits around the spin-down stationary state (marked by the blue arrow) when the gate pulse is applied, 
(d) Same as in (a) but with a fivefold increase in the intensity of the gate. The time-averaged condensate emission is almost 
circularly unpolarized when the pulse is applied {{sz) = 0). (e) Numerical calculations for (d), showing almost equal amplitude 
oscillation of the spin-up and spin-down densities, (f) Trajectory of pseudospin, showing that with a strong gate pulse large 
orbits commence around the two spin-up and spin-down stationary states, (g) Average circular polarization when the gate is on, 
at different gate intensities. When an opposite polarization gate is applied, {sz) converges to zero, whilst it remains unchanged 
for the same-spin gate. The orange line is a guide to the eye. (h) Numerical calculations of (sz) vs normalized pump intensity 
(n is the condensate occupation). 


left-circularly polarized (‘asymmetric pumping’). Once 
the condensate is created, we excite the sample reso¬ 
nantly from the back side with an oppositely circularly 
polarized (cr+) gate. FIG. 7(a) shows the intensity of 
and cr_ components of the condensate emission during 
this gating. We see a reduction in intensity for the com¬ 
ponent which is opposite to the gate laser polarization 
and an increase for the same spin polarization compo¬ 
nent. To account for this in the theory, we add a new 
term to Eq. 2, corresponding to the resonant laser: 

Gi = - l9+i{s)'ip+i - U'y- 

- i(ai|V»+i|2 +q;2|V’-iP)V’+i (7a) 

^-1 = - 5fl-i(5')V’-i - 5(7 - i£)^+i 

- 5(ai|V’-il^ + a2|^/’+iP)V’-i, 

where AIl{t,to,St) = AH{t — to)H{to St — t) uses 
the Heaviside step function H to give a square pulse 
with amplitude A which starts at time to and lasts 
for St, and ujg is the excitation frequency. To ac¬ 
count for the elliptically-polarized pumping we modify 


g{S) to g±i{S) = F — W±i + gS. Numerical calcu¬ 
lations [FIG. 7(b)] show the case W-i = 1.09 VE+i, 
A = 5 X 10 where n is the occupation of the conden¬ 
sate, and ujg = —0.3^. The dotted line marks the time 
average of the oscillations. The condensate pseudospin 
vector processes around the stationary state pseudospin 
vector at a frequency of ci;L/27r ^ 10 GHz in a limit cycle 
[see FIG. 7(c)]. In the case where the condensate is highly 
spin-polarized and the gate is at resonance {ujg = 0), the 
oscillation frequency is equal to the self-induced Larmor 
precession frequency ujl = Ig (Appendix F). This sets 
the fastest possible spin dynamics in the system, and does 
not depend on the gate intensity. Here, the condensate 
pseudospin oscillates faster than our detection time reso¬ 
lution and consequently we only see the average effect in 
FIG. 7(a). The coherent driving reported here strongly 
depends on the detuning of the gate laser frequency rel¬ 
ative to that of the condensate. The resonance width at 
which we can drive the condensate is determined to be 
10-20 iieV, as explained in Appendix E. 

In the case of EIG. 7(d) the intensity of the gate laser 
is increased by five times with respect to that in (a). 
Under these conditions, the time average of the circular 
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FIG. 8. (a) Intensity of the spin-up and -down states as a function of time. Resonant pulse G, with opposite circular 

polarization (cr+) switches the circular polarization of the condensate, (b) Despite the slow rise time of our AOM 100ns), 
the condensate switches to the opposite spin state in less than 10 ns (the limit of our detection), as soon as the intensity crosses 
the spin-switching threshold, (c) Numerical calculations for (a), showing that the spin switches with a gate pulse with intensity 
6 orders of magnitude smaller than that of the condensate. The evolution of the pseudospin is shown in (d). The two attractors, 
marked by the blue and orange arrows, are the two stationary s^i spin-polarized states. 


component of the condensate emission becomes almost 
unpolarized. If the amplitude of the resonant excita¬ 
tion is large enough, the condensate pseudospin can cross 
over to the second attractor (the spin-up state) and form 
a large trajectory which encompasses both spin-up and 
spin down stationary states with a characteristic period 
doubling [55] [FIG. 7(e,f)]. The dotted lines in (e) mark 
the time average of the oscillations. FIG. 7(g) shows the 
power dependence of the average circular polarization as 
a function of the intensity of the gate. As the oppositely 
polarized gate intensity is increased, we see the conden¬ 
sate average circular polarization converges to zero. Here, 
we have the gate resonant with the condensate. However, 
the coherent driving strongly depends on the detuning of 
the gate laser frequency to that of the condensate. The 
resonance width where we can drive the condensate is 
determined to be 10-20 peV, as explained in Appendix E. 
That theory and experiment agree regarding the observed 
average polarization over the whole range of gate powers, 
and also the observation of a resonance (see Appendix E) 
strongly suggests that our description in terms of coher¬ 
ent driving is valid. 

In summary, the fixed points of broken parity symme¬ 
try (s^,s^) become unstable due to small perturbations 
and convert to limiting circles around stationary states. 
The precession is linear for small gate amplitudes but it 
becomes nonlinear for large amplitudes. This is mani¬ 
fested by oscillation of the condensate parameters, occu¬ 
pation and polarisation with time. In this linear regime, 
the gate pulse appears to act as if it induces an effective 
magnetic field, around which the condensate polarisation 
processes. It is worth noting that this system can exhibit 
chaotic behavior with increased amplitude of the pulse. 
Increasing the amplitude first results in period doubling, 
and eventually leads to chaos (Eeigenbaum scenario) [55]. 
The ability to coherently control the spin of the conden¬ 


sate suggests its utility for computational operations. 


C. Linearly-polarized pump: spin switching 

We now explain the case when the nonresonant pump 
is linearly polarized. In this case there is no external 
‘force’ to drive the condensate to a specific spin state. As 
a result, once the circular component of the condensate 
pseudospin crosses zero, it falls into the opposite spin 
state’s attractor. EIG. 8(a) shows a realization which 
lasts for 4ps with the condensate first initialized in the 
spin down state (s^). We resonantly excite the conden¬ 
sate with an oppositely polarized gate (cr+). This causes 
the switching of the polarization of the condensate to the 
same polarization as that of the gate. The gate, which is 
1.5 ps in duration, is then turned off but the condensate 
remains in the switched polarization because the pump 
is linearly polarized and the symmetry is not explicitly 
broken. This shows that we are capable of manipulating 
the polarization of the condensate on demand, while also 
reinforcing our observation that the condensate picks a 
polarization spontaneously when it is formed. It should 
be noted that the state of the condensate does not change 
when the circular polarization of the gate is the same as 
that of the condensate. 

If the symmetry breaking was somehow set by the pa¬ 
rameters of the experiment, as in the previous example 
where the pump was slightly elliptical, the condensate 
would switch back to its original state after the gate 
beam was turned off. Instead, a very small gate field 
in the cavity can initiate coherent switching. The upper 
panel in EIG. 8(a) shows the transmitted gate intensity 
when the pump is blocked. Indeed, the gate power be¬ 
fore the cavity is 600 pW, which is 50 times weaker than 
the pump power. However, only a small fraction of this 
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resonant beam couples to the cavity in our setup. There¬ 
fore, a suitable comparison of the intensity of the gate 
to that of the condensate is to compare the transmitted 
intensity of the gate to the intensity of the condensate 
in the same direction. We find that the gate intensity is 
more than 60 times weaker than the condensate intensity, 
when the gate laser frequency is tuned to the bottom of 
the polariton dispersion (accounting for the blueshift of 
the condensate). 

By estimating the occupation number of the conden¬ 
sate (see Appendix D) we find that only 13 polaritons are 
enough to reverse the spin state of the condensate. The 
condensate switching time is less than 10 ns (our detec¬ 
tion limit), and is 10 times faster than the switching time 
of the AOM [FIG. 8(b)]. This important fact shows that 
the condensate does not follow the optical gate adiabat- 
ically. Simulations of the spin-switching with minimum 
resonant gate intensity are shown in FIG. 8(c). The po¬ 
larization of the gate is opposite to that of the condensate 
and we have assumed symmetric pumping rates for the 
spin-up and spin down condensates i.e. VF+i = W-i. 
The condensate polarization reverses within ^ 200 ps 
once the gate is applied and remains in that state after 
the pulse is turned off during continued spin evolution 
[FIG. 8(d)]. Switching requires a minimum gate inten¬ 
sity, to twist the condensate pseudospin onto the equator 
in the Poincare sphere. This sets a threshold for the 
gate power. By measuring the gate flux, and using the 
switching time of 10 ns, an upper bound of the minimum 
energy for switching the state of the condensate is found 
to be ^1 X 10“^^ J, comparable to state-of-the-art opto¬ 
electronic switches with similar speeds [56]. We empha¬ 
sise, however, that the theoretical limit for the minimum 
switching energy is 50 times smaller because the spin dy¬ 
namics of the system (^ 200 ps) is 50 times faster than 
our detection limit. This polaritonic system is then an 
extremely low-power switch. 

Multistability has indeed been demonstrated before in 
resonantly pumped condensates [11, 26]. However, there 
are several major differences between our system and 
that of the Deveaud group. In Deveaud’s experiments 
the physical process causing multistability is the nonlin¬ 
ear nonradiative losses in the polariton gas due to the 
formation of biexcitons. These losses are only significant 
when the polariton gas energy is close to the biexciton 
energy (< 2meV) [57]. In contrast here, the bistability 
that we observe is present > lOmeV below the biexci¬ 
ton energy. There is no phase transition in Paraiso et al. 
work [26] on multistability and symmetry is not spon¬ 
taneously broken. In contrast, we observe spin symme¬ 
try breaking while the power thresholds of the left- and 
right-circularly polarized states are the same, allowing us 
to observe spontaneous magnetization. In the Deveaud 
multistable system, in order to switch the polarization 
one has to inject an opposite-spin polariton density equal 
to the density difference of spin-up and spin-down po¬ 
laritons. In contrast, in our experiments, the conden¬ 
sate switches with a gate intensity 60 times weaker than 


the condensate. Theoretically, the gate intensity can be 
much weaker and the reason that we do not see switching 
at even weaker powers experimentally is due to spin fluc¬ 
tuations in the condensate, as discussed below. Finally, 
we note that other microstructured systems such as mi¬ 
crodisk lasers can also show bistability. The whispering 
gallery modes in microdisk lasers can have clockwise or 
counterclockwise propagating lasing modes [58]. When 
the coupling between the two modes is small, a cross¬ 
gain saturation causes one of the modes to dominate. 
The system can therefore operate in “flip-flop” mode and 
the two modes can be excited and switched with weak 
optical pulses [59]. 


D. Thermal noise: spin flipping 


(3) f experiment ] theory 



FIG. 9. (a) Measured spontaneous spin flip rate at dif¬ 

ferent temperatures of the sample. Above the temperature 
Tc — 15 K, the flip rate increases exponentially. The solid 
line shows the theoretical fit. (b) Calculated dependence of 
flip rate with normalized noise amplitude. The flip rate in¬ 
creases exponentially as the noise amplitude increases (light 
region), and eventually becomes linear (dark region). Com¬ 
parison with (a) shows that the measured flip rates remain 
below this linear regime. 

Above the spin bifurcation occupation threshold {Sc)^ 
any small perturbation induces a ‘restoring force’, which 
drives the condensate toward the attractor once the per¬ 
turbation is stopped. This restoring force keeps the con¬ 
densate around the attractor for small spin fluctuations. 
However, for sufficiently large perturbations, the conden¬ 
sate can flip to the other attractor with the opposite cir¬ 
cular polarization. In the experiment, for sample temper¬ 
atures near 5 K the observed spin-polarized states remain 
stable for many seconds (longer than the stability of our 
experiment can be maintained). However, increasing the 
temperature above Tc ^ I5K ± I induces spin flips in 
the condensate during a measurement time window of 
tm = 1.5 ms. The measured rate of the spin flipping 
beyond is a nonlinear function of the sample temper¬ 
ature as shown in FIG. 9(a). It is important to note 
that we do not have a real threshold for the occurrence 
of spin flipping here. In fact, the spin-flip rate increases 
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exponentially with temperature right from T = 0, un¬ 
til it becomes significant in the finite-time measurement 
window. We can account for this phenomenon by adding 
a thermally induced noise to our theory. This thermal 
noise /(t), which is similar to the Johnson noise, is in¬ 
cluded in Eq. (2) to give Langevin-type equations: 

^+1 = - y{S)'ip+i - 5(7 - 
i’-i =-Kj -is)'ip+i 

where fa{t) with a = ±1 is a realization of Gaussian 
random processes with zero mean (/a(^)) = 0 and J-like 
two-point correlation function 

{fcr{t)fcr'{t')) = 0 , = 2D5a,cT'5{t “ ^')- 

(9) 

At finite temperature T the intensity of the noise can 
be written approximately as + aT), with the 

W-contribution being a shot noise from the reservoir [16] 
and the thermal part aT defined by spin-flip polariton- 
phonon scattering. 

The flip rate vs noise amplitude using stochastic sim¬ 
ulations [FIG. 9(b)] reveals an Arrhenius-like increase at 
a critical threshold followed by a crossover to a linear 
regime. This model gives an excellent account of the dy¬ 
namics [FIG 9(a)]. For temperatures beyond 18K we 
reach the time resolution of our detection. As a result, 
we cannot completely span the crossover to the linear 
spin-flip regime in our experiment. A fit of the simula¬ 
tion results to the experimental data in FIG. 9(a) gives 
a = 0.17ps“^K“^, which sets the dependence of spin 
flip rate on D(T). We can then explore how this noise 
perturbs the spin system given by equations (8). 

Inclusion of noise produces a spin flip rate that can 
overcome the effective spin potential barrier (see Ap¬ 
pendix F). For the case in which the circular degree of po¬ 
larization is high, the spin-flip process can be considered 
as a one-dimensional Kramers transition. The spin-flip 
rate Rk can then be estimated as 



We note that for 5 = 30peV and the condensate occu¬ 
pation n = 800, the zero-temperature shot-noise spin- 
flip rate (set by W) is negligible for our observation 
timescales (10 s). The critical temperature Tc{tm) given 
by Rk{Tc) = '^/tm depends dramatically on the mea¬ 
surement window time tm- While the phonon-polariton 
interaction a in our system gives Tc(1.5 ms) = 15 K, at 
lower temperatures the condensate spin lifetime rapidly 
exceeds the stability time of our experimental apparatus 
(many seconds). Modifying the phonon-polariton inter¬ 
action a thus has an enormous effect on the spin stability. 
Finally, we note that in a similar fashion to thermal noise, 
an overlapping reservoir can also induce spin noise in the 


condensate. The spin noise causes condensate spin flips, 
which result in the reduction of the time-averaged cir¬ 
cular polarization. Theoretically, this can be studied by 
introducing a noise term similar to Eq. 9, but instead of 
depending on temperature the noise intensity depends on 
the overlapping reservoir density [60]. 


V. CONCLUDING REMARKS 

In summary, we showed how spin can emerge sponta¬ 
neously in nonresonantly-pumped polariton condensates. 
We found that for trapped condensates, in the case where 
the pump is linearly polarized, parity symmetry is spon¬ 
taneously broken by spin fluctuations at the onset of con¬ 
densation. Fluctuations are amplified by nonlinearities in 
the condensate formation due to the energy and lifetime 
splitting of the linear polarization components, produc¬ 
ing a spin-up or spin-down condensate. The symmetry 
can be explicitly broken by applying a slightly elliptically 
polarized pump, which increases the likelihood of form¬ 
ing condensates with the same spin as the pump. In the 
case where the pump is linearly polarized, we switched 
the condensate state using a 60-fold weaker resonant gate 
pulse with an opposite circular polarization. This situ¬ 
ation changes when the pump is elliptically polarized, 
where instead of switching, the condensate pseudospin 
processes around stationary states in limit cycles. Fi¬ 
nally, we showed how thermal excitations can induce spin 
flips with a rate that increases exponentially with sample 
temperature. 

We demonstrated here one way to explicitly break sym¬ 
metry utilizing elliptical polarization pumping geome¬ 
tries. One could also break the symmetry by introducing 
a magnetic field to split the energy of the spin-down and 
spin-up condensates. Alternatively, the pumping symme¬ 
try could be broken with spin current injection. These 
experiments thus exhibit rich physics with potential ap¬ 
plications in sensing. 

The observation of spontaneous discretized spin- 
polarized states also has interesting consequences in the 
physics of condensate lattices. The possibility of shared 
reservoirs, and a Josephson type tunneling [61] between 
adjacent sites, could provide new phenomena previously 
unobserved in driven bosonic systems. Magnetic phase 
transitions, geometric frustration and spontaneous pat¬ 
tern formation of spin in lattices, domain formation, 
topological spin insulators, and topological defects are a 
few examples of magnetic systems that could be studied, 
all within a highly-controlled bosonic many-body system. 

While the condensed polariton lattice resembles nano¬ 
magnet arrays [62, 63], it has the inherent advantages of 
tunable nonlinearity, longer spin relaxation time, ps re¬ 
sponse, rapid optical addressing and manipulation, and 
adaptable scalability. The spin-polarized state at zero 
magnetic field is retained for many seconds, 10 orders of 
magnitude longer than the condensation time, making it 
a suitable candidate for optical spin-based memories. 
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Spin switching with only a fraction of the condensate 
density, which is a direct result of the nonlinearity in our 
system, can be used for low power optical and electrical 
sensing and spin switches. Finally, the possibility of co¬ 
herent driving allows the realization of superpositions of 
spin-up and spin-down states, which are the key require¬ 
ments for quantum information processing. 
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Appendix A: Experimental methods 

The cavity’s top (bottom) distributed Bragg reflector 
(DBR) is made of 32 (35) pairs of Alo. 15 Gao. 85 As/AlAs 
layers of 57.2 nm/65.4 nm. Pour sets of three 10 nm GaAs 
quantum wells (QW) separated by 10 nm thick layers of 
Alo. 3 Gao. 7 As are placed at the maxima of the cavity light 
field. The 5A/2 (583 nm) cavity is made of Alo. 3 Gao. 7 As. 
The microcavity sample is chemically etched from the 
substrate side to form 300 pm diameter membranes al¬ 
lowing optical access from the back of the sample for res¬ 
onant excitation [EIG. 1 (a)]. The sample shows conden¬ 
sation under nonresonant excitation [64]. The excitation 
laser is a single-mode GW TiiSapphire, which is ampli¬ 
tude modulated using an AOM with a rise time of 100 ns. 
To pattern the pump intensity a spatial light modulator 
was used [48]. 


Appendix B: Strain-induced linear polarization 
splitting 

The measured energy of the ground state X- and 
V-polarized photoluminescence far below threshold is 
shown in EIG. 10 (a). The energy splitting of the lin¬ 
early polarized modes varies across the sample surface, 
reaching its maximum at the edges and minimum at the 
center of the membranes. This observation suggests that 
the splitting is correlated to the level of strain across the 


microcavity structure, since the latter is expected to pos¬ 
sess a similar spatial dependence: namely to be strongest 
at the boundary between etched and non-etched regions 
and relax towards the central parts of the membranes. 
Note that we observe all the same phenomena even in 
unetched samples, as strain is universally present from 
the III-V heterostructure growth. 


Strain-induced splitting of the initially degenerate po- 
lariton states at /c = 0 into orthogonal linearly polarized 
modes has been demonstrated in previous works [45, 46]. 
Both the excitonic and the photonic parts of the polariton 
can be affected by strain to produce such an anisotropy. 
In the former case, the initial splitting of the bright ex- 
citon states due to exchange interactions are enhanced 
by strain-induced mixing of the heavy and light hole 
valence bands, thereby reducing the symmetry of the 
QW [45, 65, 66 ]. In the latter case, strain induces a small 
birefringence in the cavity and/or DBRs, hence lifting the 
degeneracy between the [110] and [110] axes [46]. 


Because of the finite curvature of the cavity stopband, 
splitting the lower polariton into two orthogonal modes 
necessarily induces a difference of their linewidths as well. 
Eor the microcavity structures studied in this work the 
polariton modes are located on the low-energy side of the 
stopband [EIG. 10(b)]. In this case the mode possessing 
higher energy will exhibit a narrower linewidth, since it is 
located closer to the stopband center, where the DBR re¬ 
flectivity is at its maximum. Correspondingly, the lower 
energy mode will exhibit a larger linewidth, since it is lo¬ 
cated closer to the edge of the stopband, where the DBR 
reflectivity starts to drop. Eor an energy splitting of the 
order of 30peV, transfer matrix calculations predict a 
linewidth difference of approximately 3peV [EIG. 10 (c)]. 


This small difference in linewidth is not resolvable with 
our instruments. However, it should be noted that any 
nonzero linewidth difference, as long as it is negative sign 
relative to the energy splitting, eventually leads to the bi¬ 
furcation of the circular polarization at a critical thresh¬ 
old [EIG. 11 (a)]. Moreover, if the ratio of the linewidth 
difference to energy splitting is kept constant (as we have 
7 = O.le according to our transfer matrix simulations), 
the condensate becomes linearly polarized at high energy 
splittings [EIG. 11 (b)]. This is indeed what we observe 
at the edges of the membrane, where the energy splitting 
is as high as 100 peV. 
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FIG. 10. (a) Energy-resolved emission of the ground state polariton for horizontal (H) and vertical (V) polarizations, split in 
energy by ~ 30peV (b) Transfer matrix calculations (red) and experimental reflectivity of sample. The lower polariton mode 
(marked by an arrow), is offset by lOmeV below the center of the stopband, (c) Transfer matrix calculations showing the 
linewidth difference versus the energy splitting. The linewidth difference is 10% of the energy splitting. 
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FIG. 11. (a) Absolute degree of polarization vs 7 for a fixed 
energy splitting oi e — 0.045 ps“^. (b) Absolute degree of 

polarization vs energy splitting £, for the case when 7 = O.le. 


matrix form read 

1 1 
— =-[w - r^s - - ri)Szaz] - -(7 - ^£)crx^ 

~ 9 [(<^1 + 0^2)5' + («! — a . 2 ) Sz ( yz \ 

(C 2 ) 

and the equations for pseudospin components are 

S^ = {w- r]S)S^ - 76' - aS,Sy, (C3a) 

Sy = {w- r]S)Sy + eS, + (C3b) 

S, = {w- iiS)S,-eSy. (C3c) 

1. Linearly polarized condensates (parity 
conserved) 


Appendix C: Variable cross spin saturation 


For the case when the same-spin and cross-spin gain- 
saturation nonlinearities are different we have: 

V >+1 [w - ^{n\ip+i\^ + {2r] - /i)|V’-iP)] V’+I 

- 1(7 - i£)V’-i - 5 + a2\'ip-i\‘^] 'ip+i, 

(Cla) 

Gi =5 [w - + i2r] - ii)\^p+i\‘^)] 

- |(7 -i£:)t/+i - i [ai|t/)_i|2 + q;2|^/’+iP] V’-i, 

(Clb) 


These particular solutions are given by 5'^ = 5';^ = 0, 
Sx = and S = (re T 7)/^- X-polarized condensates 
(upper sign) can exist for VF > F + 7. V-polarized states 
(lower sign) appear for VF > F — 7. Considering their 
stability with respect to small fluctuations: Sy = y, Sz = 
z, Sx = ^So + and S = So T x, with 5'o = ('^ T 
we have linearized equations 

X = —tjSqx^ (C4a) 

y = ± 7 ^ + (5 ± aSo)z, (C4b) 

i = (±7 - (/i - y)So)z - ey. (C4c) 

Taking x^y^z (x we see that fluctuations in x always 
decay, while (C4b,c) produce the equation for the Lya¬ 
punov exponent A, 

(A T if + (m - ri)So{\ T 7 ) + ± ctS) = 0. (C5) 


where w = W—V and we have now two saturation param¬ 
eters, T] and /i. The saturation is controlled by individual 
occupations of circular components when y = 2r], and 
the saturation is controlled by the total occupation when 
/i = 77 . In general, y ^ y ^ 2r]. Equations above in the 


The stability of the V-polarized condensate is lost 
when one root of this equation crosses zero. This corre¬ 
sponds to the critical occupation Sc and critical pumping 


^c = 


7 ^ + 5^ 

[ae - (m - v)l] ’ 


= r - 7 + TjSc. (C 6 ) 
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The stability of X polarized state depends on the in¬ 
terrelation between 7 and — For 11 = r] this state 

is always unstable. However, if /i > 77 the stable X po¬ 
larized state can be formed for large enough condensate 
occupations (far above the threshold). 


optical transmission efficiencies, and R = 1.9 x 10^^ s“^ 
is the spatially-integrated count rate of the CCD. Insert¬ 
ing these values in Eq. (Dl) gives the particle number 
n ^ 800, at P = l.TPth* 


2. Elliptically polarized condensates (parity 
broken) 


Appendix E: Resonance width 


These solutions with ^ 0 are given by 




= [5^ + (re - r]S){w - /rP)] , 

ae 

{w- r]S) 


= ±e\ 


S^-Sl 


-h (re — iiSy ' 


Substitution into (C3a) gives 


(C7a) 

(C7b) 

(C7c) 


(/i - r])s [e^ + (re - r]S){w - /rP)] 

+ a‘^s{w — fiS)S + 0^7 [5^ -h (re — /rP)^] = 0. 


(C8) 


The positive root of this equation for S should be taken. 
Also, it is necessary to satisfy the condition \Sx\ ^ S', 
which gives re ^ {W 2 — 7) = Wc- This means that the 
weak lasing solutions appear continuously from the Y- 
polarized solution at the critical pumping IF2. 

The stability of weak lasing states can also be lost. 
Numerical analysis shows the following typical scenario 
of evolution of the condensate polarization state with in¬ 
creasing re for /r > 77. First, the Y linearly polarized 
state is formed. Then it transforms into a weak lasing, 
elliptically polarized state. The stability of the latter is 
also lost with increasing re, resulting in some irregular, 
quasi-chaotic dynamics and/or in the oscillatory motion 
of the pseudospin vector. Finally, for large re the stable 
X polarized state is formed. 


Appendix D: Estimate for particle number 

The condensate particle number is experimentally 
measured by: 


4>r 


(Dl) 


where 4> is the photon flux, r = 10 ps is the polariton 
lifetime, and ICp = 0.4 is the photon Hopfield coefficient. 
The photon flux is measured by: 


aR 

v 


(D2) 




gate detuning (geV) 


FIG. 12. (a) The reduction of the spin-down intensity versus 

the detuning of the gate laser frequency is shown. A sharp 
resonance is visible when the laser is in resonance with the 
condensate, (b) The numerical calculations show a similar 
resonance effect, with the center of the peak slightly redshifted 
due to the linear polarization energy splitting £. 

We study the reduction of the opposite component 
of the condensate circular polarization to the gate laser 
[marked by As- in FIG. 7(d)] as the ‘detuning’ of the 
gate varies, while the gate power remains constant. In 
order to change the detuning of the gate laser frequency 
with respect to the condensate, instead of changing the 
frequency of the laser, we tune the condensate frequency 
by changing the pump intensity. The trapped condensate 
‘blueshifts’ as the pump intensity is increased, mainly due 
to the repulsive interaction between the polaritons. This 
contrasts with the case of untrapped condensates, where 
interactions of polaritons with the excitons in the reser¬ 
voir is the source of blueshifts [32, 33, 67]. In the 4 spot 
trapped geometry, we have a blueshift of 6 peV/mW and 
as a result we can tune the condensate energy accurately 
with respect to that of the gate. FIG. 12(a) shows As- 
with respect to the detuning of the gate laser. We observe 
a sharp resonance, with a full width at half maximum 
(FWHM) of 11 ± 6peV. FIG. 12(b) shows the theoreti¬ 
cal curve versus the gate detuning {hujg), which resembles 
a Lorentzian profile with a linewidth of 17 peV. Theoret¬ 
ically, the linewidth of the resonance is defined by the 
noise in Eqs. (8). The FWHM of the Lorentzian reso¬ 
nance peak is D/{1 — sl)n. We note that the resonance 
frequency ujc is slightly redshifted, due to the effect of the 
linear polarization splitting and the decay rate splitting. 


where a = 3.5 e“/count is the photoelectron sensitivity 
of the CCD, 77 = 0.0021 is the total detection efficiency 
including the camera quantum efficiency and the total 


LUc 


1 

2 


(o^i + 0^2)5^ ~ 


{jSyS.^eSSx) 

(52 + SI) 


(El) 
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Appendix F: The Kramers flip rate 


U{S,) = 50 (1^2 - 52 \nS,). (F3b) 


From Eq. (8) we can obtain the spin vector equations: 

= -g{S)S, -jS- aS,Sy + (Fla) 

4 = -g{S)Sy + 5 ^, + + Fy{t), (Fib) 

S, = -g{S)S,-eSyFF,{t), (Flc) 


This expression for the effective potential U{Sz) is not 
valid for small Sz^ where it diverges logarithmically. The 
top of the barrier should be cut off when Sz be comes 
comparable to | {Sy) |, that is for Sz — Sq^J go /e. The 
value of poten tial o n the top of the barrier is then Uh — 
—goS^ \n{So^/go/e). The bottom of the well is positioned 
dit Uo = goSQ{^ — In So)- As a result, the spin should 
overcome the barrier 


where the correlators of real-function noise Fi(t) are 
{Fi{t)Fj{t')) =2DS6ij6{t-t'), ij = x,y,z. (F2) 


Here we consider limiting case when two parity break¬ 
ing states are formed near the north and south poles of 
the Poincare sphere, i.e., \Sz\ ^ \Sx,y\- The spin com¬ 
ponents of fixed states in this limit are Sxo — 

Syo — Sz — ^So, where So is the root of 

S = ^e/ag{S). In what follows we also denote go = g{So) 
and by assumption go 5 , 7 . Being excited away from 
the fixed state, the spin exhibits fast self-induced Larmor 
precession and slow relaxation. The precession frequency 
is cj = o:Sz: so that u — ^ejgo near the stationary states. 

The spin should be driven by noise into the equatorial 
plane {Sz = 0) in order to flip. This Kramers prob¬ 
lem can be simplified if we perform averaging over fast 
precession of the spin vector. This treatment is valid 
as long as \Sz\ ^ \Sx,y\‘ Consider the motion in the 
north hemisphere. We assume that the number of po- 
laritons does not change during the flip, i.e., S is fixed 
to So^ and we consider the case of large occupations, 
In^o ^ 1. Omitting the ^-terms from Eqs. (Fla,b) we 
find that for a given value of 5';^ > 0 the averages over 
one cycle of the other two components are {Sx) = —eja 
and (Sy) = —^So/cxSz- Then, from (Flc) we obtain the 
equation for slow evolution of Sz 

^ = - 9 oS.+go^^+F,it) = (F3a) 


AU = Ub-Uo = -goSi 


In ( - 1 - 1 

go 




(F4) 


Using the known result for the Kramers first passage time 
in the one-dimensional problem [68], we obtain the spin- 
flip rate 


Rk = ^ exp 
Zir 


AU 


Ro 


50^0. In 


2D 


go 


(F5) 


^0 = ^. 

o^go 


This expression assumes goSo/D ^ 1 and £ ^ go- 
The pre-exponent Ro cannot be written exactly by this 
method, since we do not know the shape of the effective 
potential near the top of the barrier. It can be estimated 
as Ro ^ 


Appendix G: Simulations and numerical parameters 

The parameters used for all the simulations are as fol¬ 
lows: 

OD simulations: 77 = 0.01 ps“^; T = 0.1 ps“^; he = 
30peV; 7 = O.I 5 ; hai = lOpeV; 0^2 = —0.5^1 

2D simulations: hai = 3peVpm^; 0^2 = —O.So^i; = 

46peVpm^; gp = gp/F, A = 0.1; m* = 5.1 x 10“^me; 
Jr = 10ps“^pm^; he = 7peV; 7 = 0.2^. 
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